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Abstract 



We study double integral representations of Christoffel-Darboux kernels associated with 
two examples of Hermite-type matrix orthogonal polynomials. We show that the Fredholm 
determinants connected with these kernels are related through the Its-Izergin-Korepin-Slavnov 
(IIKS) theory with a certain Riemann-Hilbert problem. Using this Riemann-Hilbert problem 
we obtain a Lax pair whose compatibility conditions lead to a non-commutative version of the 
Painleve IV differential equation for each family. 



1 Introduction 

Let us denote with (p n ) n eN the classical Hermite polynomials such that deg(p„) = n and 



As we know from the pioneering work of Gaudin and Mehta [16], the so-called Hermite kernel 



describes the statistical properties of the eigenvalues of a random matrix in the space of (n x n) 
Hermitian matrices equipped with the measure fi(M) := e _Tr ( M ^dM, where dM denotes the 
standard Haar measure. More precisely, /i(M) induces a measure on the space of configurations 
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K n (x,y) := ^2pk(x)p k (y)e 
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fc=0 



of n points on the real line, hence a determinantal point process whose particles are given by the 
eigenvalues of M, and whose correlation functions Pk{xi, ■ ■ ■ ,Xk) are given by the formula 

Pk{xi, ...,x k )= det(K n (xi,Xj))i J=1 . 

In particular, the last particle distribution F(s), describing the probability that the largest eigen- 
value is smaller than s, is given by the Frcdholm determinant F(s) = dct(Id — Xs^n), where we 
denoted with Xs the indicator function of the semi infinite interval [s, oo), and K„ is the integral 
operator whose kernel is K n (x,y) (for a very nice introduction about determinantal random point 
processes and random matrices see [15]). A remarkable connection between F(s) and the Painlcvc 
IV equation has been discovered in the nineties by Tracy and Widom. Namely, in [21], they proved 
that the log-derivative R(s) := d s log(F(s)) solves the sigma-form of Painleve IV equation 

(R") 2 + 4(R') 2 (R' + 2n) - A{sRl - R) 2 = 0. 

The aim of this article is to extend this result to the case of Christoffel-Darboux kernels associated 
to Hermite-type matrix-valued orthogonal polynomials (MOP). To that purpose we will show first 
how to obtain double integral representations of some examples of MOP. 

First, we need to introduce some notations. In what follows we consider N fixed and we use 
preferably boldface letters to denote matrices, and standard font for scalars. We also use In 
for the (N x TV) identity matrix, omitting the explicit reference to its dimension when it cannot 
lead anyone into confusion. Let C be any piecewise smooth oriented curve. A weight matrix 
W = (Wij)fj =1 : C — > GL(N,R) on the curve C is a positive definite matrix at any point of the 
curve with finite moments. We say that a matrix function F belongs to the space L 2 v (C,R NxN ) if 

J F(z)W(z)F T (z)dz < oo, 

where the superscript "T" denotes, as usual, the transpose 1 . In the above definition we mean that 
the integral is finite entry by entry. In the case when the weight matrix W is the identity matrix 
Jjv we will just write L 2 (C, R NxN ) . This induces a matrix-valued inner product for any two matrix 
functions F, G e L 2 w (C,R NxN ), denoted by 

(F,G) W = J F(z)W(z)Q T (z)dz. (1.1) 

This is not an inner product in the common sense, but it has properties similar to the usual scalar 

1 /2 

inner products. It is also possible to define a scalar norm of a matrix function F by Tr ((F, F)w) 
(see [7]). Therefore L 2 v (C,M. NxN ) with this norm is a Hilbert space and (1.1) is the inner product. 

A sequence {P n )nef-i 01 orthonormal MOP with respect to a weight matrix W is a sequence of 
matrix polynomials satisfying 

degP„ = n, (P n ,P m ) w = I N S„ mi V n, m e N. 

We will mainly work with MOP on the real line, therefore C = R. Work in the last few years 
has revealed a number of explicit families of MOP on the real line. In many cases they are joint 

1 More generally one could consider matrix- valued complex weights, and in this case you have to substitute the 
transpose with the Hcrmitian conjugate. 
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eigenfunctions of some fixed differential operator with matrix coefficients independent of the degree 
n of the polynomials. This study was initiated in [8], but nontrivial examples had to wait until 
[9, 11]. These examples are the matrix analogue of the classical families of Hcrmite, Laguerre and 
Jacobi polynomials. 

Given a complete orthonormal family of MOP (P n )neN m ^vv(^' ^ NxN ) the Christoffcl-Darboux 
(CD) kernel is defined as (see, for instance, (2.26) of [12]) 

71-1 

K n (x,y) :=Y, P k(y)Pk(x), x,yeR. (1.2) 

k=0 

We observe immediately the following properties: 

1. K n (x,y) = K^(y,x). 

2. VFe4 (M,M ArxAr ) , F{y) = (F(x),K n (x, y)) w (reproducing kernel property). 

3. K n (x,z) = {K%(z,y),K n (y,x)) w . 

Observe, in particular, from the second equation, that we are thinking about the kernel as an 
integral operator acting on the left 2 for functions in L 2 ^ (R,M WxAr ). 

In Section 2 we find double integral representations of the CD kernel (1.2) (or rather a slight 
modification of (1.2)) for two Hcrmite-type families of MOP, already introduced in [9]. We first 
find integral representations of the families (P n )neN using the corresponding second-order differ- 
ential equation that they satisfy. To the best of our knowledge, this is the first time that integral 
representations of families of MOP are studied in detail. 

Secondly, in Section 3, we study the Fredholm determinant det(Id — x s K„) of the integral 
operator K n with kernel K n {x,y). Namely we identify, using the theory of integrable operators a 
la Its-Izergin-Korepin-Slavnov [14], this Fredholm determinant to the isomonodromic tau function 
associated to a specific Riemann-Hilbert problem (sec Appendix A and B). This Riemann-Hilbcrt 
problem, through a standard procedure, can be reduced to one with constant jumps, leading to a 
certain Lax system of equations. The compatibility conditions of the Lax system give a couple of 
matrix ordinary differential equations. Combining these two equations we obtain a non-commutative 
version of the derived nonlinear Painleve IV differential equation for each family, see Theorems 3.3 
and 3.4. The contents of this section are very close, in spirit, to the one obtained by one of the 
authors and Marco Bertola in [6]. 

Finally, in Section 3.1.1, we give a symmetric formulation of the non-commutative Painleve IV 
equation. The non-commutative Painleve II equation used in [6] has been introduced by Retakh 
and Rubtsov in [18], where the authors obtained this equation as a reduction of a non-commutative 
analogue of Toda equations. Hence it appears desirable to verify if also the non-commutative 
Painleve IV equations obtained in this article can be written, analogously, as reductions of some 
suitable non-commutative Toda-type equations. The content in Section 3.1.1 goes in this direction. 

2 It is possible to work with a CD kernel acting on the right, but in that case we have to consider a different 
inner product defined by (F,G) W := J R G T (x)W(x)F(x)dx. Now the CD kernel will be defined by K^(x,y) := 
^Z 1 P?(x)P k (y) = K n (y,x). 
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2 Hermite-type MOP and related CD kernels 

In this section we consider a couple of examples of MOP already introduced in [9]. Both are 
orthogonal with respect to Hermite-type weight matrices of the form 

W(x) := e- x2 T(x)T T (x), x E R, 

where T is certain matrix polynomial. Consider an orthonormal family of MOP (P n )neN- Now 
define, for every n e N, the orthonormal function 

*„(*) :=e-*'' 2 P n (x)T(x). 

(<S>n)neN is a family of matrix functions orthonormal with respect to the identity matrix, i.e. 

($ n ,* m ) /jv = I N S nm , Vn.meN, 

where the inner product (•, •) is defined by (1.1). In our examples the family (<J>„) ne N will always 
be complete in the space L 2 (R,M. NxN ^ (see Section 6 of [13] for more details). 
The Hermite-type CD kernel is then defined, slightly modifying (1.2), as 3 

71-1 

K n {x,y):=Y,*k{v)*k{x). (2.1) 

fc=0 

In order to find double integral representations of these two examples of matrix-valued ker- 
nels (2.1) we will use the fact that the corresponding MOP are eigenfunctions of a second-order 
differential equation of the form 

P:(x)F 2 (x) + P' n {x)F 1 {x) + P n (x)F = T n P n (x), (2.2) 

where F% , F\ and Fq arc matrix polynomials (which do not depend on n) of degrees less than or 
equal to 2, 1 and 0, respectively, and the eigenvalue T n is a symmetric matrix. 

2.1 The first example 

Let A N be the (N x N) nilpotent (A% = 0) matrix 

N 

A N :=J2^E t , i+ i, CjgR, (2.3) 
»=i 

where Eij = (Si r S S j)^ s= i is the elementary matrix with 1 at entry (i, j) and elsewhere, and Jn 
the diagonal matrix 

JV 

J N :=Y,(N-i)E^. (2.4) 

i=l 

3 The kernel defined in this section, indeed, is equal to the one defined in the introduction up to a conjugation by 
e _ TT(x). In order to keep the notation simple we used, nevertheless, the same symbol. 
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Again, we will sometimes remove the dependence of N and write A — An and J = Jn, whenever 
there is no confusion about the dimension of the matrices. A and J satisfy the algebraic relation 
[A,J] = -A. 

Let W be the following weight matrix 

W(x) =c- x2 e Ax c ATx , xel, (2.5) 

already introduced in [9]. Observe that e Ax is an upper triangular matrix polynomial of degree 
N -1 (since A N =0). 
The family of MOP 

P n (x) = e- A2 / 4 P„(x), (2.6) 

where (P„)„gN denotes the monic orthogonal family with respect to (2.5), satisfies a second-order 
differential equation as in (2.2) (see, for instance, Section 4 of [13]) where 

F 2 (x)=I, F 1 (x) = -2xI + 2A, F (x) = A 2 - 2J, r„ = —2nl - 2 J. 

The family (2.6) is not orthonormal, but it will be normalized later for the computation of the CD 
kernel (2.1). In the following, given a matrix M, we will use the standard notation z M := c Mlogz , 
and the branch cut of log z is chosen to be the real negative axis. 

Theorem 2.1 Let (P n )neN be the family of MOP defined by (2.6). Then there exist suitable 
constant matrices C n and D n such that 

P n (x)c Ax = jf z - J C n z J c- z2+2zx ^, (2.7) 
\{x)e Ax = e x ' J w J D n w' J e w2 - 2xw w n dw, (2.8) 



and 

Pn 

Ix 

where the contour 7 encloses the origin, closes at —00 and it is traversed in a counterclockwise 
direction while X := L + iR, and L > is chosen so to have no intersection between 7 and I (see 
Figure 1). 

Proof: For the first integral representation (2.7) we observe (using the formula e~ Ax Je Ax = J+xA) 
that the functions Y n {x) — P n {x)e Ax satisfy the following differential equation 

Y;'(x) - 2xY^(x) - 2Y n (x)J + 2{nl + J)Y n (x) = 0. (2.9) 

Let us find solutions of (2.9) of the form 



l^*) = / V n (z)e 2zx dz, 



where r] is some contour in the z-plane. Substituting this expression into the differential equation 
(2.9) and integrating by parts, it is easy to see that (2.9) holds if the two following conditions are 
satisfied 

V; L (z) = -{2z+ ^±A) V n (z) + - z [V n (z)J - JV n {z)] , (2.10) 
zV n (z)e 2zx = 0. (2.11) 
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General solutions of homogeneous first-order matrix differ- 
ential equations of the form 

X'{z) = A(z)X(z) + X(z)B(z) 

are easy to obtain considering solutions X(z) = Xi(z)CX 2 (z), 
where C is any constant matrix, such that X[(z) = 
A(z)X 1 (z) and X' 2 {z) = X 2 (z)B(z) (see for instance [22], 
Lemma 30.1). Applying this to the equation (2.10), there ex- 
ist, for every n > 0, constant matrices C„ such that the general 
solution is of the form V n (z) = Vn,i{z)C n V n ^{z), where 

=-(2zI+^I+±>)V nil (z), VUz) = V n , 2 (z) J -. 



Therefore 



V n (z) 



C n z - 



■n+l 




Figure 1: Contours 7 and I. 



gives a solution for (2.10). Now we observe that, choosing 
77 = 7 (the same contour as in the case of scalar Hermite 
polynomials), the expression on the left hand side of (2.11) is 

nothing but the residue at infinity of the expression zV n (z)e 2xz , which is clearly zero because of 

2 

the term e~ z . Hence 

dz 



Y n (x) := <p z J C n z J e 



yn+1 



is a solution of (2.9). In particular, since P n (x)e Ax satisfies a differential equation as in (2.9), 
there will be suitable constant matrices C n , which depend on the family (P n ) ne tq, such that we get 
the integral representation (2.7). 

For the second integral representation (2.8) we follow a similar argument, just observing that 
the functions Z n (x) — e~ x P n (x)e Ax satisfy the differential equation 



Z;'(a;) + 2xZ' n (x) - 2Z n (x)J + 2((n + 1)/ + J)Z n {x) = 0, 
and considering Fourier type solutions of (2.12) of the form 

/oo 
v n {ty 
-00 



(2.12) 



.,2ixl. 



dt. 



Then (as before) we find some first-order matrix differential equation for V n (z) and we conclude 
that 

V n (t) = t J C n t- J e- t2 t n , 



where C n ,n > 0, are certain constant matrices. Hence 



Z n (x) 



t J C n t~ J e 



-t -\-2ixtj.n 



t n dt. 



G 



In particular, since e~ x P n (x)e Ax satisfies a differential equation as in (2.12), by the change of 
variables t — > iw, followed by a shift iw — > iw + L, we conclude that there will be suitable constant 
matrices D n = i n+1 i J C n i~ J , depending on the family (P n ) ne ^, such that (2.8) holds. Q.E.D 

Remark 2.1 Observe that the constant matrices C n and D n obtained initially from the differential 
equations (2.9) and (2.12), respectively, have many degrees of freedom. It is exactly the choice 
of the family {P n ) n &i which allows us to have a unique representation of C n and D n , as we will 
see below. The meaning of this is because there are many families of matrix polynomials satisfying 
(2.9) or (2.12), but not all of them are orthogonal with respect to some weight matrix. 

Corollary 2.1 The Hermite-type CD kernel (2.1) associated with the MOP with weight (2.5) can 
be written as 

K n {x,y) =c (z 2 -2/ 2 )/2 f dw I dzM n (z,w)e w2 ~ 2xw ~ z2+2zy , (2.13) 
Jx J 7 

where 



M n (z,w) := z- 



C k'<- 



fell w 



W 



Lfe=o 



i 



-w 



(2.14) 



and \\Pk\\w denotes the matrix-valued norm with respect to (1.1). 

Proof: It is just enough to use the formula 

* n {x)= C - x2 l 2 \\P n \\^P n {x)e Ax 

and then the integral representations (2.7) (with the substitution x — > y) and (2.8) in the definition 
of the matrix-valued kernel (2.1). Q.E.D 

The coefficients C n and D n are not easy to obtain for our example since there are no general 
structural formulas like the norms or the coefficients of our family of MOP for any general size N. 
We will study in detail the case N = 2 in the next subsection. 



2.1.1 A detailed study of the case N = 2 

Let us consider the (2 x 2) case related to MOP with respect to the weight 



W(x) = e~ x e A2X e A * x , A 2 := 



v 




v e 



We have many structural formulas for this example (see [10]). In particular, if we denote 



7 2 

In 



1 + -v z 



we have that the polynomials (P„)„ 6 n in (2.6) have diagonal norms 



1 2 _ nl V^ ( ll+l 



n\\W — 



2 n 







i/ri 



The coefficients C n and D n of the integral representations (2.7) and (2.8) are given by 



1 v(n+l)/2 



(2.15) 



(2.16) 



7 



and 



1 

u/tt 



1 V 

— 1/7 2 



(2.17) 



These coefficients are computed directly from the integral representations (2.7) and (2.8), since 
we have that the family (P„)„ e N can be given in terms of scalar Hermite polynomials (see [10]). 
Therefore (2.16) is obtained from (2.7) using Cauchy's residue theorem and (2.17) is obtained 
from (2.8) using standard manipulations of the Fourier transform of Hermite polynomials. 

Proposition 2.1 The Hermite-type CD kernel (2.13), for N = 2, can be written as 



K n {x,y) = tt^U^ 2 -^/ 2 { dwldz z J *B n z- J 'w J *B- 1 w- J * 



— 2xw — z +2zy+nlog(w/z) 



W — Z 



where 



B n := 



1 

nv 



, det(B„) = 7*. 



(2.18) 
(2.19) 



Proof: Computing (2.14), using (2.15), (2.16) and (2.17), we have that 



M n (z,w) 



+ 



(7)' 



w — z w — z 



zjll - 1) + w 

,2 



wi n 



'ST., 



V 



(7)' 



2 7 2 ZW 



+ 



(7)' 



w — z w — z 



will -l) + z 



z ln 



J 



Then, using Cauchy's residue theorem, we observe that the diagonal terms —l/(w — z) and the — v 
in the upper right corner of M n (z, w) do not give any contribution. Hence we have 



K n {x,y) = 



(2tu) 2 



-y 2 ^ 2 [ dwl 

JX J -1 



dz M n (z, w)- 



— 2xw — z +2zy+n\og(w/z) 



W — Z 



where 



/ z(il-l) + w 



M n (z,w) = 



wX 



V 2 7 2 



(z — w) 



-z(w-z) 

In 

will - 1) + Z 

Zln 



and finally it is easy to see that 

M n (z,w) = z J *B n z- j2 w J *B- 1 w- J \ 
where B n is given by (2.19). Q.E.D 



8 



Remark 2.2 Also for the case N = 3, the Hermite-type CD kernel (2.13) has a representation 
analogue to (2.18), for a given (3 x 3) constant matrix B n (we do not report here the computations, 
since they are completely analogues but more cumbersome). We conjecture that for any N > 1 and 
a given n > 0, there always exists a constant matrix B n such that 

9 2 , N f f e w 2 -2xw-z 2 +2zy+nlog(w/z) 

K n {x,y)=- -c (x ~y )/ 2 / dw idz z jN B n z- ,N w jN B- 1 W - jN .(2.20) 

(2m) 2 Jx J~f ui-z 

2.2 The second example 

Let W be the following weight matrix 

W{x) =e- x2 e Bx2 e BTx2 , (2.21) 

where B = A(I + A)" 1 and A is the (N x N) nilpotent matrix (2.3). The family of MOP 

P n (x) = [(I + A)- 1 / 2 } 2n+1 P n (x), neN, (2.22) 

where (P„)„ s n denotes the monic orthogonal family with respect to (2.21), satisfies a second-order 
differential equation as in (2.2) (see Section 5 of [13]) with 

F 2 (x)=I, F 1 (x) = 2x(2B-I), F (x) = 2(B-2J), T n = -2nI-4J, 

where J is the diagonal matrix (2.4). Also in this case, in complete analogy with the first ex- 
ample, we can deduce two different integrable representations of the polynomials (P n )neN and, 
consequently, an integral representation of the related Hermite-type CD kernel. 



Theorem 2.2 Let (P n )neN be the family of MOP defined by (2.22). Then there exist suitable 
constant matrices C n and D n such that 

P n (x)e Bx2 = j z- 2J C„z 2J e- 2 + 2 --g T , (2.23) 

and 

P n (x)c Bx2 =e x2 J w 2J D n w- 2J c w2 - 2xw w n dw, (2.24) 

where the contour 7 encloses the origin, closes at —00 and it is traversed in a counterclockwise 
direction while 1 := L + iR, and L > is chosen so to have no intersection between 7 and X (see 
again Figure 1). 

Proof: This theorem is proven exactly in the same way as Theorem 2.1, using as starting point the 
formula e' Bx2 Je Bx2 = J + x 2 (B - B 2 ) instead of e~ Ax Jc Ax = J + xA. Q.E.D 

Corollary 2.2 The Hermite-type CD kernel (2.1) can be written as 

K n {x,y) = c ( x2 ~v 2 y 2 J dw j> dzM n (z,w)e w2 ' 2xw - z2+2zy , (2.25) 



where 

"n-l 

,2J 



M n (z,w) = z A 



Y^Clz- 2J \\P k \\^w 2J D k 

k=0 



-w- 2J . (2.26) 



and \\Pk\\\v denotes the matrix-valued norm with respect to (1.1). 
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As before, the coefficients C n and D n are not easy to obtain in general for this example. We 
will focus on the case N = 2. 

2.2.1 A detailed study of the case N = 2 

Let us consider the (2 x 2) case related to matrix-valued polynomials orthogonal with respect to 
the weight 

W(x) = e-^e^V^ 2 , B 2 = A 2 . 
Again, from [10], we have many structural formulas for this example. If we denote 



n(n — 1) 



we have that the polynomials (P n )neN m (2.22) have diagonal norms 

||2 _ ( 51+2 



\\r n \\w~ 2 n y l/6 2 n 

The coefficients C n and D n of the integral representations (2.23) and (2.24) are given by 



Cn 



2 n+1 m 



and 



1 



v(n + l)(n + 2) 

1 V 

n{n-l)v . . 



4# 



V5l 



(2.27) 



(2.28) 



(2.29) 



(2.30) 



which can be computed in a similar way as the first example. 

Proposition 2.2 The Hermite-type CD kernel (2.25), for N = 2, can be written as 



K n {x,y) = jT—^e^ >/ 2 dwidz z 2J -B n z- J *w J *B n w- 2J - 
where B n is a (2 x 3) matrix and B n a (3 x 2) matrix given by 



—2xw—z +2zy+n\og(w / z) 



W — Z 



(2.31) 



Br, :- 



5l + i 



vn(n + 1) nv 



B„ :- 



\ 45l +1 25 2 l+1 S 2 n J 

B n is a right inverse of B n , i.e. B n B n = I 2 . 



vn(n — 1) 1 
V 451 I J 



(2.32) 
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Proof: Computing (2.26), using (2.28), (2.29) and (2.30), we have that 



/ 



M n {z,w) 



{2m) 



;\2 



(!)-'-' , (!)-', (i) 



W — Z 



w6 2 



v —w — z + 



^ +i hi) , "("-^(i) \ y - 1 y y 

4(5:;, ,z 2 w ■ 4<5 2 zw 2 w - z z^ 2 ,, z<5 2 



rl T c-9 



tf 2 



V 



V 



^ 2 n+1 z 2 w 



Then, using Cauchy's residue theorem, we observe that the diagonal terms —l/(w — z) and the 
— w — z in the upper right corner of M n (z, w) do not give any contribution. Now taking out the term 
(™) (w— z)^ 1 of the above expression we observe that it can be written as z 2j2 B n z~ j3 w j3 B n w~ 2j2 
where B n and B n are given by (2.32). The equality B n B n = I2 is a consequence of the definition 
of S 2 t , see (2.27). Q.E.D 

Remark 2.3 As before, we also analyzed the case N = 3. The Hermite-type CD kernel has a 
similar representation (2.31) for certain constant matrices B n and B n of size (3 x 5) and (5 x 3), 
respectively, satisfying B n B n = J3. We conjecture that for any N > 1 and a given n > 0, there 
always exist constant matrices B n and B n of size (N x (2JV — 1)) and ((27V— 1) x N), respectively, 
satisfying B n B n = In and such that 



K n {x,y)=- 



2(wi) 



- dw &>i 

Jx J~i 



dz z 2Jn B n z~ j2N - l w j2N - 1 B nW - 2jN> - 



— 2xw~z +2zy+nlog(w/z) 



W — Z 



-. (2.33) 



3 Hermite-type kernels and non-commutative Painleve IV 



Let us consider a (TV x N) matrix-valued kernel of the general form 



K n (x,y) = T^e^ 2 )/ 2 f dwcfdz B n (z)B n (w)- 

(27TZ) Z Jx J 



-2xw~ z 2 +2zy+n \og(w / z) 



W — Z 



(3.1) 



where 3 n and 3 n are, respectively, an [N x p) and a(px N) square-integrable matrix function 
such that & n (z)B n (z) — In- We think of K n {x,y) as a kernel of an integral operator K n acting 
on the left for matrix functions, namely 



[K n F}(x):= f F(y)K n (x,y)dy VFeL 2 (E,l 



,NxN 



)• 



The Hermite-type CD kernels studied in Section 2 are particular cases of this one. 

In the following Xs will always denotes the indicator function of the interval [s, 00). We start 
with the following theorem: 



Theorem 3.1 Consider, given sel, the operator '. 
with kernel (acting on the right) 



: L 2 ( 7 UI,: 



oNxN 



)^i 2 (7Ul, 



pNxN 



K n . s (w, z) :-- 



>log( Z ) 



2-ki(w — z) 



c -Tr+2 2 +™l°g«S( z )B( w ) 
2iti(w — z) 



X 1 {w)xi{z). 
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The following equality between Fredholm determinants holds: 

det(Id - x s Kn) = dct(Id - J 

Proof: Using the isomorphism L 2 ("/ U I, M. NxN ) = L 2 (-f,l 
in matrix form as 



n,s ) 
NxN\ ^ t2( 



L 2 (l,R NxN ) we can write K„„ 








F s " 




G 






where the operators F s : L 2 (j,R NxN ) — -> L 2 (l,R NxN ) and G : L 2 (l,R Nx ") — ► L 2 (-f,R NxN ) 
are defined respectively by the kernels 

T s {\z) := 



e -^-+2s( Z -A)-nlog(z) 



27ri(A- z) 



Xi(A)x 7 (2;), 



Q{zM := X.WxxH- 

Let us also introduce the Hilbert-Schmidt operator s written in matrix form as 



(3.2) 
(3.3) 



K 






s • 


id 


F s " 





Id 






-F s ' 









Id 


-F s " 




' Id — F s o G 





-G 


Id 




-G 


Id 



Through the identity 

(id-f; iS )(id-f„, 5 ) = 

we get the following chain of equalities (see the Appendix B for the definition of det 2 ): 

det(Id - f „,„) = det 2 (Id - f „,,) = dct 2 (Id - K^ a )det 2 (Id - f „ >s ) = 

= det 2 (Id - ¥ s o G)c- Tr ( F = oG ) = dct(Id -F.oG). 

Now, using the formulas (3.2) and (3.3), we deduce that F s o G : L 2 (l,R NxN ) — > L 2 (l,R NxN ) 
has kernel explicitly given by the convolution 

dz c w2 - z2+2 < z -^ +nio ^- nio ^B n {z)B n {w) 
{2m) 2 ~ (X-z)(w-z) 

Finally we conjugate (F s o G) with the Fourier transform T such that 

(TF)(x)= ( 4^c- 2Aa; F(A), {T- 1 G){X)= ( -^Le 2Xx G(x). 

Then we obtain (using Cauchy's residue theorem) that the kernel C(x,y) associated to (ToF s oGo 
T -1 ) is equal to 

,,2 2, 



(3.4) 



(3.5) 



C(x,y) 



dX 

— ( 



e 2A(x- S ) 1 dz / dyJ 



e w -z +28z-2wy+nlog(w)-nU>g(z)fg n ( z jjg n ( w } 



(2ni) 2 (X- z)(w- z) 
if x < s 



, 2 .x 2 I dw (£ dz B n (z)B n (w) 



—2xy—z -\-2zx-\-n \og(w/z) 



W — Z 



(3.6) 



if x > s. 
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and this latter, up to a conjugation with the operator of multiplication by e ^ , is equal to Kj(x, y). 
Hence (3.6), together with (3.5) and (3.4), gives (3.1). Q.E.D 

We now introduce a Riemann-Hilbcrt problem related to the kernel K n ^ s through the IIKS 
theory. In the following we denote 

6 n (X,s) := A 2 - 2As + nlog(A). 



Problem 3.1 Find the sectionally analytic function T(A) G GL(N+p, C) onC\{-fi)l} such that 

r+(A)=r_(A)(J-G(A)), XenUj, 



r(A) = + S + A ^°°- 



with 



G(A) := 



c e "( A ^)B„(A) 




Xx(A) 





_ e -MA, s ) B T (A) q 



X 7 (A)- 



Theorem 3.2 The Fredholm determinant det(Id — x s M. n ) is equal to the isomonodromic tau func- 
tion (B.3) related to the Riemann-Hilbert problem 3.1. Hence, in particular, we have that 

d\ 



5 s logdet(Id- Xs K„)= / Tr(n 1 (A)(a A r_)(A)H(A)) 



2m 



where we denoted 



E(A) := 8 S (I - G(X))(I - G(A))- 1 = -d s G(X)(I + G(A)). 



(3.7) 
(3.8) 



Proof: 



Let's define the two ((N + p) x N) matrices written block-wise as 

e- 2sA /jvXi(A) 



/(A) := 



9(A) := 



1 

2^i 



-e-V B T (A)x7(A) 
e-^+ 25A -" lo s( A )JjvX 7 (A) 



c V+„log(A ); g n(A)xi(A) 

It is straightforward to verify that the following two equations are satisfied: 



-, G(X) = 2mf(X)g'(X). 



A — /i 

Hence, using Theorem B.l together with Theorem 3.1, we conclude that, denoting by tjmu the 
Jimbo-Miwa-Ueno tau function related to the Riemann-Hilbert problem 3.1, we have 

TjMu(s) = det(Id - K n . s ) = det(Id - x s K„), 

and, as a consequence, equation (3.7). The second equality in the equation (3.8) comes from the 
fact that g T (X)f(X) = 0. Q.E.D 
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3.1 Non-commutative PIV 



We now specialize to the particular case in which N — 2 and K„ is the Hermite-type CD kernel 
(3.1) associated with the MOP studied in Sections 2.1.1 and 2.2.1. 

In the first case we have B n := z j2 B n z~ j2 and B n (w) := (£?„(u;)) _1 , where B n is the constant 
matrix (2.19). In the second case we have B n (z) := z 2j2 B n z~ Ja and B n (w) := w j3 B n w~ 2j2 , 
where B n and B n arc defined in (2.32). Let us introduce, for these two different (sub)cases, the 
two matrices 



T A (\) :- 



T B (A) := 



MM 

2 



I 2 - J 2 log(A) 








e n (x,s) 



^i/ 2 -2J 2 log(A) 




2 I 2 - J 2 log(A) 




and 



J 3 log(A) 



In the following we use the curly brackets to denote the anti-commutator between two matrices, 
i.e. {x, y} :— xy + yx. Moreover we denote with a prime the derivative with respect to s. 

Theorem 3.3 Let B n be as in (2.19) and T{\) be the solution of the Riemann-Hilbert problem 
3. 1 with jump 



G{\) 



e e ^ x ^\- j2 B~ T X j2 




Xi(A) + 











_ e -e n (x,e) x -j 2B T X j 2 



The 



d s logdet(Id-x 5 Kn) = Tr((r 1 ) 22 - (I^u), 
where K n is the integral operator with kernel K n (x,y) given by (2.18). 
Moreover *(A) := r(A)e TA ^ A ^ satisfies the Lax equations 
, , ; h - J2) A" 1 



X 7 (A). (3.9) 



(3.10) 



(A -)/, + ((= - 

2y~ 1 z + (y-V - y-'uz)^ 1 -(A - s)I 2 - (|/ 2 + J 2 - y^zy) A" 1 



9,* = 



-AI 2 y 
-2y- x z XI 2 



where 



* ■■= -(ri)'n, 

y := -2(ri)i 2) 

k tx:-(r 1 )' 12 (r 1 )r 2 1 + 2 S / 2 . 



(3.11) 
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The compatibility conditions give the following coupled system of ODEs: 

u' = -u 2 + 2su + 4z - 2nl 2 + V A , 



2u'z + 2uz' - 2sz', 



(3.12) 



where V A ■= 2[J 2 ,y]y 1 . 



Combining these two equations we obtain a non-commutative version of the derived PIV equation, 
in the form 

u"'+ [u", u) - 4(n + 1 + s 2 )u'- 2 ({«', u 2 } + uu'u) 

(3.13) 

+6s{u', u} + Au{u - s) + (V A - 2(uV A )Y + 2sV A = 0. 

In the second case we will be dealing with rectangular matrices. Hence in the following, given 
an (N x p) rectangular matrix M, with N < p and linearly independent rows, we denote with Aft 
the right inverse 4 of M defined as M f := M T (MM T y 1 . 

Theorem 3.4 Let £?„, B n be as in (2.32) and T(\) be the solution of the Riemann-Hilbert problem 
3. 1 with jump 

c^^A-^BjA^ 

' Xx(A) + 



G(A) 













. e -9„(A, S ) A -2J 2 _ B T A J 3 



X 7 (A). 



Then 



a s logdet(Id- Xs K„) =Tr(^(r 1 ) 22 - (ri)nj, 
where K„ is the integral operator with kernel K n {x,y) given by (2.31). 

Moreover \I/(A) := r(A)e Ts ' A ^ satisfies the Lax equations 



(3.14) 



<9a* = -4.* 



(A - a) h + 



- z) Io- 2 Jo I A 



2ytz + (yV _ yt uz ) A -i -(A - s)I 3 - (|/ 3 + J 3 - y 1 zy) A" 1 

-AJ 2 y 
-2yU XI 3 



where 



z:=-(ri)' n) 
y:=-2(ri)i 2) 

[ t t: =(r 1 )' 12 (r 1 ) t 12 + 2 S / 2 . 



4 In fact, the definition of a right inverse of an (N X p) rectangular matrix M is not unique. For any invertiblc 
(p X p) matrix C, a right inverse of M can be defined as Aft := CM T (MCM T )- 1 . Therefore wc will have 
eventually a family of equations of the form (3.15). However we normalize all the computations assuming C = I. 
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The compatibility conditions give the following coupled system of ODEs: 

u' = -u 2 + 2su + Az -2nI 2 + V B , 



2u'z + 2uz' - 2sz', 



where Vb ■= AJ 2 — 2yJ 3 y^. 



Combining these two equations we obtain a non- commutative version of the derived PIV equation, 
in the form 

u"'+ [u", u]- 4(n + 1 + s 2 )u'- 2 ({«', u 2 } + uu'u) 

(3.15) 

+6s{u', u} + Au(u -s) + (V B - 2(uV B ))' + 2sV B = 0. 

Remark 3.1 Before going into the proof let us remark that both the equations (3.13) and (3.15) 
are equations for just one variable, namely y; u and V (either Va or Vb) being functions of y. 
The reason why we claim that this equation is a non-commutative version of the derived PIV is 
that, in both cases, if we assume that all the variables commute, we get the equation 

u" - Au - 6u 2 u + Uu'u - Anu + 4w 2 - Asu - As 2 u = 0, 

and the reader can easily verify that this latter equation is the derivative of the standard PIV equation 

u" = ^+ 3 u 3 - Asu 2 + 2(s 2 + 1 + n)u - — . 
2u 2 u 

Proof of Theorem 3. 3: 

Given T(A) solution of the Riemann-Hilbcrt problem 3.1 with G{\) as in (3.9), we have 



Tr (r: 1 (A)9 A r_(A)H(A)) ^ = - res Tr (n 1 (A)9 A r_(A)9 s T A ) , (3.16) 

yUX Z7TZ A — oo 

where the (formal) residue above simply stands for minus the coefficient of the power A -1 in the 
asymptotic expansion of the argument. The formula (3.16) can be proven using Cauchy theorem 
and goes back to the article of Palmer [17]. A very precise detailed derivation is given in [4], section 
5.1. Direct application of the formula (3.16) and the equation (3.7) yields the equation (3.10). 

Next we observe that ^(A) = r(A)e TA ( A ) solves a Riemann-Hilbert problem with constant jump, 
hence both (Sa*)* -1 and (<9 S 1 4') , I , ~ 1 are meromorphic functions on C* so that, in particular, 
considering the singular behavior at at oo, we get 

(dx*)*- 1 = A=XA 1 + A + \~ x A-i 

(9,*)*- 1 = U = XUi+Uo- 

Let us start with A; in order to compute the coefficients Ai,Aq, we compute the first terms of the 
asymptotic expansion at infinity of (dx^)^^ 1 giving immediately 5 



Ai = Ii ® er 3 , Ao = -sl 2 ® er 3 + [Tx, I 2 <g> 0-3] = sh ® <?3 ' 



-2(ro 12 

2(ri) 2 i 



3 We denote with 0-3 the standard Pauli matrix 0-3 = diag(l, —1). 
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The equation for A-i is slightly more complicated: we use again the expansion at infinity and we 
get 

ft 

A-i = [r 2 , i 2 ® <t 3 ] + [i 2 ® <t 3 , ri]ri + - j 2 ® <t 3 - j 2 ® j 2 - [r l7 s / 2 ® <x 3 ]. (3.17) 

In order to simplify the expression above we use the fact that the the A _1 -tcrm of (<9 s \l r )\f'_ 1 
term is identically zero, giving 



leading eventually to 



ri = [r 2 ,j 2 ®o- 3 ] + [i 2 ®ff 3 ,ri]ri, 



A-i = r'j - [T u sl 2 ® cr 3 ] + -I 2 <g> cr 3 - J 2 



(3.18) 



The block-diagonal part of the equation above gives 

(r 1 )' 11 = 2(r 1 ) 12 (r 1 ) 21 , (r^ = -2(ri) 21 (r 1 )i 2) (3.19) 

and so, combining (3.17), (3.18) with the definitions (3.11) and using the relations (3.19) we get 



An 



-Si 2 V 

2y~ 1 z sJ 2 



n r 

-z + -i 2 - J 2 



uy 

' 2 



(^-i) 



21 



In order to get the missing term (.4._i) 2 i we use the Lax equation 

d s A-dxU = [U,A]; 

the (1,1) term giving 

(.4.-1)21 = y~V - y^ 1 uz. 

In this way we arrive to the stated expression for A. The computation for U, since we have no 
singularities at the origin, is simpler and the asymptotic at infinity of (<9 S \J') 1 I ,_1 gives immediately 



Ui = -I 2 <E)cr 3 , U 



2(ri)ia 
-2(ri) 21 



y 

-2y~ x z 



The first equation in (3.12) is just the (1, 2) (block) entry of the equation. For the second one 
we start observing that the off-diagonal terms of the relation (3.18) gives the equations 



(ri)i 2 = -2(r 2 )i2 + 2(r 1 ) 12 (r 1 ) 22 , (r^ - 2(r 2 ) 21 - 2(r 1 ) 21 (r 1 ) 11 . 



(3.20) 



Then we go on analyzing the asymptotic expansion at infinity of {d\^)^ 1 ; the A 2 -term (which 
is identically zero) gives the equation 



i\ = [r 3 , i 2 ® <r 3 ] + [ s i 2 ® <r 3 , r 2 ] 



Ti, -I 2 ® er 3 - J 2 <g> I 2 



(3.21) 



+[r 1 ,s/ 2 ® <T 3 ]ri+[ri,i2 ® <r 3 ]rf +[/ 2 ® cr 3 ,r 1 ]r 2 +[i 2 o T 3 ,r 2 ]r 1 . 
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The (1, l)-entry of this equation gives 

(ri)n = -2 S (r 1 ) 12 (r 1 ) 21 - 2(r 1 ) 12 (r 1 ) 21 (r 1 ) 11 

-2(r 1 ) 12 (r 1 ) 22 (r 1 ) 21 + 2(r 1 ) 12 (r 2 ) 21 + 2(r 2 ) 12 (r 1 ) 21 , 



(3.22) 



and combining (3.20) with (3.22) we get the second equation in (3.12) (the computation is lengthy 
but completely straightforward). The equation (3.13) is obtained simply expressing z in function 
of u using the first equation in (3.12) and then substituting in the second one. Q.E.D 

Since the proof of the Theorem 3.4 is formally identical to the one of Theorem 3.3 we will not 
write it here. One has just to substitute J 2 with J3 where appropriate, and substitute y~ l with 
y^. Indeed, going through the computations of the proof above, is easy to verify that we use just 
the property yy^ 1 = I 2 , while y~ 1 y = I2 is never used. 

3.1.1 A symmetric formulation of PIV 

While formulas (3.13), (3.15) establish a direct connection with the classical PIV equation and 
Tracy-Widom results [21], in view of equations (3.10) and (3.14) it would be desirable to write 
a system of ODEs in which both the entries (ri)n and (ri) 22 appear explicitly as dependent 
variables. To this aim the symmetric formulation of PIV given (in the scalar case) by Aratyn, 
Gomes and Zimerman in [2, 3] is particularly suitable. A similar description can be given in our 
non-commutative case. We introduce new notations for the entries of the Lax matrices so to match 
with the cited articles. 

Theorem 3.5 Let *(A) := T(A)e T ^ (A) be as in Theorem 3.3 and denote 



q 


■= 2^)12, 


r 


:= -2(ri) 2 i 


PR 


:= ^rOn, 


PL 


:= -4(ri) 22 



Then the related Lax matrices A. and U reads 

U = 



XI 2 q 
r XI 2 



(A - s)I 2 q 

-r -(X-s)h 



+ 



1 

4A 



p' R + 2n/ 2 - 4 J 2 Asq + 2q' 

Asr - 2r' ~p' L - 2n h - 4J 2 
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while Lax equations read 



p R = 2sqr + q'r - qr' , 
p L = 2srq + rq' - r'q, 



sq' + ^{~q" + 2qrq) = (1 + 2n)q + 2[q, J 2 ] 



sr ' _ _( r " _ 2rqr) = (-1 + 2n)r + 2[r, J 2 \. 



(3.23) 



Theorem 3.6 Let *(A) := r(A)e TB < A ' be as in Theorem 3.4 and denote 



q 


■= 2(ri) 12 , 


r 


:= -2(r!) 21 


PR 


:= ^rOn, 


PL 


:= -4(ri) 22 



T/ien ffte related Lax matrices A. and U reads 

U = 



A = 

while Lax equations read 



(A - s)I 2 q 
-r -(X-s)I 3 



+ 



1 

4A 



AJ 2 q 
r AJ 3 



p' R + 2n/ 2 - 8 J 2 4sq + 2q' 

4sr - 2r' -p' L - 2nl 3 - 4J 3 



p R = 2sqr + q'r - qr' , 
p L = 2srq + rq' - r'q, 



(3.24) 



-sq' + -(-q" + 2qrq) = (1 + 2n)q + 4[q, J 2 ), 

sr' - ~{r" - 2rqr) = (-1 + 2n)r + 2[r, J 3 ]. 

Remark 3.2 TTie system of equations above (3.23) and (3.24) are the non-commutative analogues 
of equations (2.4) in [2], which are equivalent to the sigma-form of PIV (see equations (2.6), (2.7) 
in [2]) 

Proof: Both theorems are proven in the same way and the proof consists, essentially, on rewriting 
the same equations as in Theorems 3.3 and 3.4 with different variables. 

In particular the third and the fourth equations in (3.23), (3.24) come from the the Lax equation 



dsA-dyU = [U,A], 
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together with the already used fact that p' R = — qr and p' L = rq (see equations (3.19)). The 
other two equations come from the (already used) equations (3.20) combined with (3.22) and 

(ri) 22 = 2 ( s (r 1 ) 21 (r 1 ) 12 + (r 1 ) 21 (r 1 ) 11 (r 1 ) 12 + (ri) 2 i(ri)i 2 (ri) 2 2 - (ri) 2 i(r 2 )i 2 - (r 2 ) 2 i(ri)i 2 ) , 

this latter coming from the (2, 2)-entry of (3.21). Q.E.D 

4 Concluding remarks 

In this paper we have shown (for the first time, to the best of our knowledge) some integral rep- 
resentations of two particular examples of Hermite-type MOP on the real line, and then deduced 
a double integral representation of the related Christoffcl-Darboux kernel K n (x,y). This was the 
starting point for the study of the Fredholm determinant det(Id — x s K„), where K n is the integral 
operator with kernel K n (x,y). Using some Riemann-Hilbert techniques, we related this Fredholm 
determinant with a non- commutative version of the Painleve IV equation. 

One interesting question that arises is how these kernels K n (x 7 y) behave if we consider scaling 
limit as n — > oo. It is very well known that, after rescaling appropriately the variables, the Hermite 
kernel converges to the Airy kernel in the following way 

lim 1 ^Hermite ( ^ + X V \ = K Ai (x,y), 

where x and y are in a bounded set (see [1] for a very nice and elementary deduction of this 
convergence, even in a more general setting than what needed here). 

In the matrix case, for the two Hermite-type kernels we study in this paper (see (2.20) and 
(2.33)) there will be a scalar behavior as n — > oo, i.e. 

lim 1 — K n [V*i+ X V2^+ 7 J— )=K Ai (x,y)I N . 

In particular, this property holds for the examples (2.18) and (2.31) (N = 2). This is easy to see 
from the integral representation (2.20) (analogous for (2.33)) using classical steepest descent meth- 
ods, as in [1]. Indeed the "matrix part" of the kernel can be written as z Jn B n (— ) Jn B^ 1 w^ Jn , 
but under the given rescaling 6 ™ — >■ 1 as n — »■ oo, therefore the matrix part converges to the identity 
matrix In as n — » oo, no matter the choice of B n . 

The Airy kernel is related, as it is well known, to the Painleve II equation [20], and this relation 
extends to the matrix case for a specific type of matrix Airy kernel [6]. For the two Hermite- 
type polynomials we study in this paper, the scaling limit of the corresponding CD kernel has a 
scalar behavior, so that the non-commutativity disappears and the corresponding matrix Painleve II 
equation is nothing but N non-interacting copies of the scalar equation. It would be interesting to 
find other Hermite-type MOP whose kernels do not behave in this way, something that in principle 
is not an easy task. This consideration, however, goes beyond the scope of this paper, and will be 
pursued elsewhere. 
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6 Rccall that the variables w and z arc rcscalcd as w = ^J~^ (l + ^173) an d z = (l + n i/7t ) > scc W- 
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A A brief reminder of regularized Fredholm Determinants 



We refer to [19] for the relevant details: we shall need only the elementary facts which we recall 
here. In general the Fredholm determinant of an operator of the form (Id — G) can be defined only 
when G is of trace class. Recall that if G is represented as an integral operator on a (separable) 
Hilbert L 2 (X, d/x) with kernel G(x,y) (we abuse notation here) then 

oo ^ „ n 

det(Id-G) = 1 + ~\ / det[G(a; i ,a; j )] iJ <„[JdM(a; i ). (A.l) 

There are other trace ideals I p , p e N, which means that G p is trace-class [19]; in particular Z 2 
consists of Hilbert-Schmidt operators. For Gel p one can define following Carleman a regularized 
determinant det p (Id — G) which has the same main property of vanishing iff the operator is not 
invertible. In particular for Hilbert-Schmidt operators one has 

dct 2 (Id - G) := 1 + V" — / det[G{xi,Xj)(l- Sij^ijKnWdfiixi), 
*=i n ' Jxn i=l 

that is, one simply omits the diagonal elements in the determinants under the integral sign. This 
determinant has the properties 

• if G is also trace-class then 

det 2 (Id - G) = dot (Id - G)c- trG . 

• if Gi,G 2 are Hilbert-Schmidt operators (and hence GiG 2 is trace class) then 

det 2 (Id - Gi)det 2 (Id - G 2 ) = dct 2 (Id - Gi - G 2 + GiG 2 ) C - tr ( GlG2 ). 

An interesting occurrence (which is used in this article) is that if G is just Hilbert-Schmidt but its 
kernel vanishes on the diagonal G(x, x) = then the series defining det 2 (Id — G) is identical to the 
regular det(Id — G). The reason for still wanting to distinguish det 2 from det in this case is simply 
that G may fail to have a trace and in a different basis the ordinary det may simply be ill-defined. 

B Integrable kernels and isomonodromic tau functions 

In this appendix we recall some basic facts about integrable kernels a la Its-Izergin-Korepin-Slavnov 
and their connections with isomonodromic tau functions, recalling in particular a theorem proved 
in [5] (see also [6]). Given a piecewise smooth oriented curve C on the complex plane (possibly 
extending to infinity) and two matrix-valued functions 

/,g:C— J-Mat pxfe (C), 

we define the kernel K as 

A — fi 
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We say that such kernel is intcgrable if f T (X)g(X) = (so that it is non-singular on the diagonal). 
We are interested in the operator K : L 2 (C,C k ) — > L 2 (C,C k ) acting on fc-vector functions via the 
formula 

(Kh)(\) = J K(\,n)h(j*)dn, 

and, in particular, we are interested in the Fredholm determinant det(Id — K) defined as in (A.l). 
The key observation is that, denoting with d the differentiation with respect to any auxiliary 
parameter on which K may depend, we obtain the formula 

aiogdetpd-A") = -Tr((Id + R)dK), (B.l) 

where R is the resolvent operator, defined as R = (Id — K)~ 1 K. Moreover R is again an integrable 
operator, i.e. 

A — fj, 

and F, G can be found solving the following RH problem: 

' r+(A) = r_(A)M(A), AeC, 

< r(A) = / + C(A- 1 ), A — >oo, (B.2) 

_ M(A) = /-27rz/(A)g T (A). 

More precisely we have the two equalities 

F(A) = r(A)/(A), G(A) = (r-f(A) ff (A). 

Now suppose, as in the formula (B.l), that the operator K (and hence the Riemann-Hilbert 
problem (B.2) depends smoothly on a certain parameter set of parameters 7 . On the space of these 
deformation parameters, we introduce the following one-form (here below d denotes a vector in the 
space of deformation parameters) 

c M (a):=^Tr(r: 1 (A)9 A r_(A)H 9 (A))^, (B.3) 
S a (A):= dM(X)M~ 1 (X) 

The definition (B.3) is posed for arbitrary jump matrices; in the case of the Riemann-Hilbert 
problem (B.2) the spontaneous question arises as to whether % in (B.3) and the Fredholm 
determinant are related. The answer is positive within a certain explicit correction term, as in the 
theorem below 

Theorem B.l ([5]) 8 Let /(A; s),g(X; s) : CxS — > Mat pX fe(C) and consider the Riemann-Hilbert 
problem with jumps as in (B.2). Given any vector field d in the space of the parameters S of the 
integrable kernel we have the equality 

Lu M {d) = aiogdct(Id -K) + H(M), 

7 In the case treated in this article, we just have one parameter, namely s. 

8 Actually the article [5] treats the case k = 1, but the proof does not change considering this more general case. 
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where lom{9) is as in (B.3) and 



H(M) := Hi(M) - H 2 (M) = J Ti{df T g + f T dg)dX - 2m J Tr(g T f'dg T f)dX. 

In the cases we treat in the article, moreover, we have H(M) = 0. Hence it is possible to define, 
up to normalization, the isomonodromic tau function tjmu '■= ex P(J w m) and this object, thanks 
to the previous theorem, will coincide with the Fredholm determinant det(Id — K). 



References 

[1] Adler, M., Ferrari, P. L. and van Moerbeke, P., Airy processes with wanderers and new 
universality classes, Ann. Probab. 38 (2010), no. 2, 714-769. 

[2] Aratyn, H., Gomes, J. F. and Zimerman, A. H., On the symmetric formulation of the Painleve 
IV equation, see arXiv: 0909.3532. 

[3] Aratyn, H., Gomes, J. F. and Zimerman, A. H., Darboux-Backlund transformations and 
tational solutions of the Painleve IV equation. Nonlinear and modern mathematical physics, 
146-153, AIP Conf. Proc, 1212, Amcr. Inst. Phys., Melville, NY, 2010. 

[4] Bertola, M., The dependence on the monodromy data of the isomonodromic tau function, 
Comm. Math. Phys. 294 (2010), no. 2, 539-579. 

[5] Bertola, M. and Cafasso, M., The transition between the gap probabilities from the Pearcey to 
the Airy process-a Riemann-Hilbert approach, Int. Math. Res. Not. 2012, no. 7, 1519-1568. 

[6] Bertola, M. and Cafasso, M., Fredholm determinants and pole-free solutions to the noncom- 
mutative Painleve II equation, Comm. Math. Phys. 309 (2012), 793-833. 

[7] Damanik, D., Pushnitski, A. and Simon, B., The analytic theory of matrix orthogonal poly- 
nomials, Surveys in Approximation Theory 4 (2008), 1-85. 

[8] Duran, A. J., Matrix inner product having a matrix symmetric second order differential op- 
erator, Rocky Mountain J. Math. 27 (1997), 585-600. 

[9] Duran, A. J. and Griinbaum, F. A., Orthogonal matrix polynomials satisfying second-order 
differential equations, Int. Math. Res. Not. 10 (2004), 461-484. 

[10] Duran, A. J. and Griinbaum, F. A., Structural formulas for orthogonal matrix polynomials 
satisfying second-order differential equations, Constr. Approx. 22 (2005), no. 2, 255-271. 

[11] Griinbaum, F. A., Pacharoni, I. and Tirao, J., Matrix valued spherical functions associated to 
the complex projective plane, J. Funct. Anal. 188 (2002), no. 2, 350-441. 

[12] Griinbaum, F. A., de la Iglesia, M. D. and Martincz-Finkelshtein, A., Properties of matrix 
orthogonal polynomials via their Riemann-Hilbert characterization, SIGMA 7 (2011), 098, 31 
pages. 



23 



[13] de la Iglesia, M. D., Some examples of matrix-valued orthogonal functions having a differential 
and an integral operator as eigenf unctions, J. Approx. Theory 163 (2011), no. 5, 663-687. 

[14] Its, A. R., Izergin, A. G., Korepin, V. E. and Slavnov, N. A., Differential equations for 
quantum correlation functions. In Proceedings of the Conference on Yang-Baxter Equations, 
Conformal Invariance and Intcgrability in Statistical Mechanics and Field Theory, volume 4, 
pages 1003-1037, 1990. 

[15] Johansson, K., Random matrices and determinantal processes, Lectures given at the summer 
school on Mathematical statistical mechanics at Ecole de Physique, Les Houches, see arXiv: 
math-ph 0510038. 

[16] Mchta, M. L., Random matrices. Third edition. Pure and Applied Mathematics (Amsterdam), 
142. Elsevier/ Academic Press, Amsterdam, 2004. 

[17] Palmer, J., Zeros of the Jimbo, Miwa, Ueno tau function, J. Math. Phys. 40 (1999), no. 12, 
6638-6681. 

[18] Retakh, V. and Rubtsov, V., Noncommutative Toda Chains, Hankel Quasideterminants and 
Painleve II Equation, J. Phys. A: Math. Theor. 43 (2010) no. 50, 505204, 13 pp. 

[19] Simon, B., Trace ideals and their applications. Second edition. Mathematical Surveys and 
Monographs, 120. American Mathematical Society, Providence, RI, 2005. 

[20] Tracy, C. A. and Widom, H., Level-spacing distributions and the Airy kernel. Comm. Math. 
Phys. 159 (1994), no. 1, 151-174. 

[21] Tracy, C. A. and Widom, H., Fredholm Determinants, Differential Equations and Matrix 
Models, Comm. Math. Phys. 163 (1994), no. 1, 33-72. 

[22] Wasow, W., Asymptotic expansions for ordinary differential equations, New York, Dover, 
1987. 



24 



